Direct detection constraints on superheavy dark matter 
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The dark matter in the Universe might be composed of superheavy particles (mass <;10 10 GeV). 
These particles can be detected via nuclear recoils produced in elastic scatterings from nuclei. We 
estimate the observable rate of strongly interacting supermassive particles (simpzillas) in direct dark 
matter search experiments. The simpzilla energy loss in the Earth and in the experimental shields 
is taken into account. The most natural scenarios for simpzillas are ruled out based on recent 
EDELWEISS and CDMS results. The dark matter can be composed of superheavy particles only if 
these interact weakly with normal matter or if their mass is above 10 GeV. 
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The dark matter in the Universe might be composed 
of supermassive particles (mass JilO 10 GeV). Although 
the leading dark matter candidate is a weakly interacting 
massive particle (WIMP), recent models for nonthermal 
production of particles in the Early Universe have broad- 
ened the dark matter mass and cross section parameter 
space. Supermassive particles avoid the unitarity mass 
bound by not being thermal relics of the Big Bang. 
There are many ways to produce supermassive particles 
in the Early Universe. The most general way is gravita- 
tional production at the end of inflation as a result of the 
expansion of the Universe 0, 0, • In this scenario the 
average particle density is independent of the interaction 
strength with normal matter. 

A reasonable assumption for the mass of such parti- 
cles is the inflaton mass scale, ~ 10 12 GeV in chaotic 
inflation models. These particles might be composed of 
exotic quarks or other new particles 0, ||J . Their produc- 
tion allows them to be electrically charged but there are 
arguments that they have to be neutral 0|. 

If the non-luminous matter is indeed composed of su- 
permassive particles its interaction strength with normal 
matter can range from weak (wimpzillas 7]) to strong 
(simpzillas .8]). Although not many experiments have 
been built to search for strongly interacting dark mat- 
ter, observations from several experiments constrain their 
mass and cross section parameter space 0, . Masses 
up to about 10 4 GeV are ruled out and several allowed 
regions were found above this mass for a cross section 
range of about 10 -32 to 10~ 15 cm 2 . Stronger constraints 
can be found in . 

Here we consider the direct detection of simpzillas. Di- 
rect detection experiments measure the energy deposited 
by a nuclear recoil produced in the particle scaterring 
from a nucleus. Their main goal is to detect WIMPs. 
We show that direct detection experiments are also able 
to probe simpzillas as a dark matter candidate. Com- 



parison of our estimated simpzilla detection rate with the 
latest CDMS [13 and EDELWEISS [HE! results, rules 
out the most natural scenarios of the simpzilla parameter 
space. 

A simpzilla arriving at an underground detector will 
have interacted many times in the Earth and, depend- 
ing on its interaction cross section with ordinary matter, 
will interact many times in the detector itself. To deter- 
mine if simpzillas can be directly detected, their energy 
degradation and range in the Earth and the experimental 
shield has to be taken into account. 

We assume that the local dark matter is composed en- 
tirely of simpzillas, with a density of 0.3 GeV/cm 3 . The 
total energy deposited (Q) in a detector is given by the 
nuclear recoil energy (En) times the number of simpzilla 
interactions in the detector Ni ~ rtffa x ffl, where is 
the detector atomic number density, a X N is the simpzilla 
nucleus cross section and I is the detector thickness. We 
assume that the simpzilla nucleus cross section is inde- 
pendent of the nuclear spin and relates to the simpzilla 
nucleon cross section (a xp ) by <j x n — a xv (Am r /m rp ) , 
where A is the nucleus atomic number, m r is the nucleus 
- simpzilla reduced mass given by m x mjq / (m x + mjv), 
m x is the simpzilla mass, m^r is the nucleus mass and m rp 
is the nucleon - simpzilla reduced mass. As m x » mjv, 
the reduced masses are simply the nucleus and the nu- 
cleon mass, respectively. 

The nuclear recoil energy is given by 

E R = J^-i = m N v 2 (1 - cos 9) (1) 

2m,N 

where q is the momentum transferred, v is the dark mat- 
ter velocity and 8 is the scattering angle in the center-of- 
momentum frame. For a given material and scattering 
angle, the recoil energy depends only on the simpzilla 
velocity. 

A simpzilla with a cross section of 10~ 26 cm 2 will in- 
teract about 10 4 times in a 1 cm thick Ge detector, while 
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it will interact only once if the cross section is 10 _30 cm 2 . 
The total mean energy depositions will be about 450 MeV 
and 45 keV respectively. 

We assume that the dark matter has a Maxwellian ve- 
locity distribution (/(«)) with an average velocity vq of 
220 km/s and include the motion of the Earth (t>@) as 
described in [l5lll6| : 
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The distance (/) traveled through the Earth to the 
detector is a function of the arrival angle between the 
simpzilla arrival direction at the Earth and the normal 
direction from the detector to the Earth's surface. We as- 
sume simpzillas are isotropically distributed in the galac- 
tic halo. The average simpzilla energy at the detector is 
given by: 



E = E c^[- 2pNA ^ N l 
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where Eq is the simpzilla energy at the Earth, p is the 
Earth density and Na is Avogadro's constant. We use 
the Earth density profile given by the Preliminary Earth 
Model and its composition found in |igj . 

Knowing the energy loss in the Earth we determine 
the maximum distance a simpzilla can travel before be- 
ing stopped. This distance is related to the maximum 
arrival angle and defines an acceptance cone. Any sim- 
pzilla which reaches the Earth at an arrival angle larger 
than this cone angle will have lost all its energy before 
reaching the detector. This cone depends on the m x and 
on a xp . The fraction of the Earth contained in the ac- 
ceptance cone (/®) will therefore define the fraction of 
incident dark matter particles that reach the detector. 

Within an acceptance cone the average simpzilla veloc- 
ity (vod) and energy (Eod) at the detector are determined. 
The fraction of these events (Jd) with energy below a cer- 
tain threshold is given by the ratio of the lower velocity 
events over the total number of events: 
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where E rnax is the maximum energy which can be de- 
tected in a given detector and E esc is the maximum 
simpzilla velocity which we assume to be 650 km/s, the 
galactic halo escape velocity. 

The deposited energy spectrum dR/dQ in the detector 
is given by: 
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where k is a parameter in the velocity distribution 
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dRi/dQ is the differential rate including the Earth's mo- 
tion 161, 
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Ro includes the fraction of the dark matter which is de- 
tectable. Besides containing the terms and fjj which 
are related to strongly interacting particles, and taking 
the energy loss into account (which makes vod < vq), 
the differential energy spectrum given by Equation [5] is 
the same as for WIMPs yjj. For a simpzilla - nucleus 
interaction the wavelength A = fi/q is smaller than the 
nuclear radius. The drop in the effective cross section 
with increasing q is described by a form factor, F 2 ( Q). 
We assume that F 2 (Q) is well approximated [lJJUJi by 
the Woods-Saxon form factor [20|. 

We consider the CDMS and the EDELWEISS exper- 
iments, which employ Ge and, in the case of CDMS, 
Si detectors operated at mK temperatures [l2l l2ll | . 
Currently, CDMS is located at a shallow depth of about 
16mwe at the Stanford Underground Facility (SUF), its 
final destination is the Soudan mine in Minnesota, at a 
depth of about 2080 mwe. EDELWEISS is located in the 
Laboratoire Souterrain de Modane at 4500 mwe. 

A particle interacting in a cryogenic Ge or Si detec- 
tor will create phonons and electron-hole pairs. While 
phonons are a measure of the total recoil energy, only a 
fraction of this energy is dissipated into ionization. The 
simultaneous measurement of the phonon and ionization 
signals results in an excellent discrimination efficiency 
against electron recoils, which are caused by the domi- 
nant radioactive background. This discrimination is pos- 
sible because the ratio of ionization to recoil energy is 
lower for nuclear recoils than for electron recoils. 

We estimate the simpzilla elastic scattering rates at 
both shallow and deep sites. We take into account the 
simpzilla energy loss in the experimental shields, made 
of Cu and Pb. For CDMS, we also estimate the en- 
ergy deposition in the active muon shield, made of plastic 
scintillator. To determine the simpzilla detection rate we 
integrate the differential energy spectrum given by Equa- 
tion EJ 

Figure^shows the detectable simpzilla rate versus sim- 
pzilla mass for various cross sections, for CDMS at SUF. 
Two effects are responsible for the kinks in these curves. 
First, the fact the experiment is at a shallow depth makes 
the distance the simpzilla travels in the Earth very small 
for arrival angles smaller than 90° and much larger above 
90°. Second, the simpzilla energy loss rate increases with 
depth, due to the increase in the Earth density. 
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The simpzilla rate at shallow sites is slightly higher 
than at deeper sites. As an example, the rates of 
10 12 GeV simpzillas arriving at a CDMS detector at 
SUF for a xp of 1(T 30 and 1(T 26 cm 2 are 27/day/kg 
and 21 x 10 4 /day/kg respectively while 4/day/kg and 
17 x 10 4 /day/kg respectively at the EDELWEISS site. 

To obtain limits on simpzilla masses and cross sections, 
we compare the predicted simpzilla detection rates with 
the background rates of the CDMS and EDELWEISS 
experiments in different energy regions. We require that 
the detectable rate translates into at least one particle 
going through the detector in the total exposure time. 
We emphasize here that our results will be conservative, 
for we are not attempting a simpzilla-specific analysis of 
the data. 

We estimate the simpzilla rate in the 10-100 keV recoil 
energy region for CDMS, and in the 30-100 keV region for 
EDELWEISS. E max in EquationHis set to 100 keV. We 
require a minimum nuclear recoil of 1 keV per interaction, 
since the ionization efficiency in Ge decreases rapidly |22^ 
below this energy. This recoil energy range probes the 
region where lower than about 10 cm . The 

a xp — 10~ 30 cm 2 curve in Figure^is representative of the 
predicted rate in this recoil energy region. For CDMS, 
we also estimate the energy deposited in the 4.1 cm thick 
plastic scintillator surrounding the detectors. Simpzillas 
with a cross section below 10 -24 cm 2 , will deposit less 
than 2MeV in this active muon shield and thus belong 
to the category of muon-anticoincident events. 

CDMS observed a total of 27 nuclear recoil events (sin- 
gle and multiple-scatters) in the 10-100 keV region for an 
exposure of 15.8kg-day (l2||. This yields a background of 
1.7 events/kg/day. Comparing the above number with 
the number of expected simpzillas in this energy region 
excludes the region labeled 'CDMS' in Figure at 90% 
CL, for cross sections below 10 -28 cm 2 . 

EDELWEISS has reported their result for WIMP 
searches based on an exposure of 11.7kg-day. It is a 
combined result from two measurements. The first [23j |. 
using a 320 g Ge detector, had an effective exposure of 
4.53kg-day. No nuclear recoils were found in the 30- 
200 keV energy range. The second measurement with 
a new 320 g Ge detector 0, had an effective exposure 
of 7.4kg-day and observed no nuclear recoil events from 
20-100 keV. Their expected background from neutron 
scatters is about 0.03 events/kg/day above 30keV [23| . 
The combined result corresponds to a rate inferior to 
0.2 kg" 1 day" 1 at 90% confidence level. Comparing our 
estimated rate for the EDELWEISS detector with the 
0.2 events/kg/day 90% CL limit for the 30-100 keV recoil 
energy region, we exclude the region labeled 'ED EL' in 
Figure |3 again for a xp below 10~ 28 cm 2 . 

For CDMS, we estimate the simpzilla rates in the 3- 
10 MeV ionization energy region, and compare it to the 
background in the same region shown in Figure 33 in 
|l2| . The muon-anticoincident background in this region 
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FIG. 1: Estimated simpzilla rate versus mass for simpzilla 
nucleon cross sections of 10~ 30 , 1(T 26 and 1(T 22 cm (from 
left to right) in the CDMS experiment at the Stanford Under- 
ground Facility. Kinks of each curve are due to the shallow 
depth and to the Earth density profile. 

is very low, due to the fact that the background from nat- 
ural radioactivity ends around 2615 keV. We approximate 
the ionization energy to one third of the recoil energy. 
The 1 keV minimum nuclear recoil energy requirement 
allows to probe o xp of 10 -26 cm 2 and lower. This com- 
parison excludes, at 90% CL, the a xp region from about 
10" 29 to 10~ 26 cm 2 labeled "CDMS" in Figure El 

Simpzilla with cross sections greater than 10~ 26 cm 2 
will interact more that 10 4 times in a 1 cm thick Ge de- 
tector, and thus deposit an energy above 10 MeV for 
a 1 keV threshold per interaction. In spite of their 
high energy, such events would nonetheless trigger the 
CDMS and EDELWEISS detectors, giving rise to satu- 
rated pulses. To obtain a conservative exclusion region 
for these cross sections we consider the total background 
trigger rate of 0.4 Hz for CDMS and the trig ger rate 
of 20 events/kg/day above an energy of 1.5 MeV [lj] for 
EDELWEISS. Since CDMS triggers on the phonon sig- 
nal |l2(, we can relax the threshold per interaction to 
1 eV. We assume a low-energy threshold of 5 keV. This 
excludes the region above 10~26 cm 2 in Figure|21 

Figure [21 also shows the predicted region to be probed 
by CDMS at the Soudan site. We assume a total expo- 
sure of 230 kg day and a muon-anticoincident nuclear re- 
coil background of 0.02 events/kg/day in the 10-100 keV 
energy region. We also assume an overall trigger rate 
which is 100 times lower than the one at SUF. 

We have investigated the direct detection of strongly 
interacting supermassive particles. We have estimated 
the differential and total direct detection rates for the 
CDMS and EDELWEISS experiments. We find that al- 
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FIG. 2: Excluded regions at 90% CL in the simpzilla mass 
versus cross section parameter space. The region labeled "ex- 
cluded" was excluded in the analysis given in [Iol| ; the regions 
labeled "EDEL" (filled) and "CDMS" (hatched) are excluded 
based on EDELWEISS [HIH and CDMS results [13 . The 
area labeled "Soudan" (within solid line) is the predicted re- 
gion to be probed by CDMS at Soudan. As a comparison, 
the region above the thick straight line shows the sensitivity 
of future, cubic-kilometer sized neutrino telescopes |24j . 



though the energy loss as well as the depletion in the 
number of simpzillas reaching an underground detector 
are substantial, the predicted nuclear recoil rates are 
much higher than for supersymmetric WIMPs. 

Comparison of our predicted rates for CDMS at SUF 
and for EDELWEISS with their most recent results 
UH EI EH rules out the most natural simpzilla scenar- 
ios. The most natural scenarios are the ones for which 
the simpzilla mass is comparable to the inflaton mass in 
chaotic inflation models (~ 10 12 GeV) and the simpzilla 
- nucleon cross section is comparable to the nucleon - nu- 
cleon strong interaction cross section. The simpzilla mass 
versus cross section parameter space which is probed by 
these two experiments is shown in Figure [2] 

The region to be tested by cubic kilometer neutrino 
telescopes such as IceCube |2_4j is also shown. These tele- 
scopes can search for secondary neutrinos produced in 
simpzilla annihilation in the Sun [8J, l25| . Although most 
of this region is excluded by direct detection experiments, 
neutrino telescopes would provide an independent confir- 
mation. 

In conclusion, the dark matter in the universe may be 
composed of superheavy particles only if these particles 
interact weakly with normal matter or if their mass is 
above 10 15 GeV or higher, depending on the simpzilla - 
nucleon interaction cross section. 
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